Co— propagating Bose-Einstein Condensates and electromagnetic radiation: Emission 

of mutually localized structures 
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Using a semi-classical model to describe the interaction between coherent electromagnetic radi- 
ation and a Bose-Einstein condensate in the limit of zero temperature, including the back action 
of the atoms on the radiation, we have analyzed the phenomenon of emission of solitary-like wave 
packets which can accompany the formation of mutually localized atom-laser structures. 
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I. INTRODUCTION 



With the realization of Bose-Einstein condensates 
(BEC) and of coherent atomic beams all the questions 
inherent to the manipulation of such systems have ac- 
quired a certain importance. In particular, early studies 
such as [l|, Q have started an interest in the manipula- 
tion of atomic structures via their interactions with co- 
herent electromagnetic radiation. Not only these studies 
could be of importance for applications such as atom in- 
tcrferometry, they also offer a possible test of the analo- 
gies between optics and quantum matter waves. In fact, 
BECs provide us with a quantum system where matter 
waves can be realized on macroscopic scales and which, 
under the approximations of zero temperature, low den- 
sities, weak interactions and within the limits of validity 
of a mean field theory, is amenable to a mathematical 
description based on the Gross-Pitacvskii equation, com- 
pletely akin to the basic equation of nonlinear optics, i.e. 
the nonlinear Schrddinger equation, The Kerr- like 
nonlinearity is given for the atoms by the atom-atom in- 
teractions. It has been demonstrated that it is possible 
to reproduce typical optical and nonlinear optical phe- 
nomena with a a BEC, from the generation of solitons 
to four-wave mixing, from parametric amplification to 
second harmonic generation, to mention only a few of 
them (for a review, see [|| and references therein). It 
is possible to push the analogy even further and con- 
sider the electromagnetic radiation as the medium that 
allows for nonlinear interactions between atoms as dis- 
cussed initially in [l[. This corresponds exactly to the 
optics case where the medium through which radiation 
propagates can bring about nonlinear effects for the elec- 
tromagnetic field: Nonlinear effects in the dynamical evo- 
lution of the atoms and of the radiation are then a con- 
sequence of the atom-light interactions. Models of this 
interaction were presented by several authors, @, @ and 
Krutitsky et al, @, gave a full derivation of the equa- 
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tions describing it starting from first principles within 
the framework of quantum field theory. This last work 
showed the emergence of a resonant nonlinear term in 
the system dynamics as a consequence of the laser-atom 
dipole-dipole interaction, besides the well known Kerr- 
like nonlinearity. The same equation for the atoms and 
consequently a coupled system of equations for the laser- 
atom system were rederived in Q , this time starting from 
a semi-classical theory. It was found there that the re- 
sponse of the " medium" , that is of the laser radiation, 
to the dynamics of the condensate could play an impor- 
tant role in the coupled evolution and even allow for the 
formation of mutually localizd atom-laser structures ca- 
pable of propagating with no changes in the atom den- 
sity and laser intensity, in spite of the assumed repulsive 
atom-atom interaction. Such solitary-like structures are 
of interest bacause of their properties of self-localization 
and robust propagation and the effects of these interac- 
tions can be seen also in relation to the creation of meta- 
lenses and comoving potentials to refocus atom waves, 
Q. Their emerging as a result of the coupling during 
propagation of atoms and laser was studied in [7| while 
the stationary states of the coupled system and their sta- 
bility properties were introduced in However, little 
is known about their actual mechanism of formation. In 
the optics case of a focusing nonlineariry, we would ex- 
pect an initial bell-shaped structure to shed away the 
radiation which cannot be accomodated and to adjust 
asymptotically to a soliton wave. It is interesting to see 
whether the same happens in the present case of coupled 
atom-laser propagation and how, since the coupling may 
lead to novel effects. We have therefore studied numeri- 
cally the process through which such coupled soliton-likc 
objects are formed evidentiating the occurrence of a phe- 
nomenon reminiscent of soliton emission in nonlinear op- 
tics, fioj . In fact, the equations analyzed here, predict 
the formation of solitary-like structures for both atoms 
and light which can move away from the region where 
they were generated. Although in our case there is no 
external trapping but only the self-consistent interaction 
of atoms and laser, these results are suggestively simi- 
lar to escaping solitons described and observed in com- 
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pletely different environments, for instance in nematic 
liquid crystals, [TT| . 

We will breafly review the basic physics of the semi- 
classical model and the limitations to be considered in 
SecHIl Sec |nl] presents an investigation of the initial evo- 
lution of the coupled system which will than be studied 
numerically in Scc lIVI 



II. SEMI-CLASSICAL MODEL AND SET UP OF 
THE PROBLEM 

The basic physics of atom-laser interactions in the sim- 
plest dipole approximation is given by photons exciting 
atoms which in turn re-emit photons absorbed by other 
atoms, th us g iving rise to a long-range interatomic in- 
teraction, [12| . Details of the semi-classical derivation of 



both the atom and the laser equations are given in 0, @ , 
here we will only briefly review the two model equations 
and re-introduce the notation. In a semi-classical deriva- 
tion, the force exerted by the light on the atoms is written 
as the gradient of a potential and this potential is used 
as the atom-laser interaction term in the Hamiltonian 
for the atoms Schrddingcr equation. Such force term is 
a generalization of the ponderomotive force and takes 
into account the possibility of an inhomogenoeus gas. 
The existence of stationary solutions is physically cru- 
cial so we will study a far-off resonant monochromatic 
field E(r, t) = Re[£(r) cxp(— iujLt)] (where £(r) is the 
complex amplitude of the laser field) . The time averaged 
force (over laser cycles) is F = y^V [|E| 2 4|-] = -VV d . 
Here e(u), n) is the medium dielectric constant with atom 
density n and is given by e(w, n) = 1 + yf^jf^;, where, 

as derived from quantum theory, a{uS) = —d 2 /h,A is the 
atomic polarizability at the laser frequency ul, with A = 
w l — w a being the detuning from the nearest atomic res- 
onance frequency u ai and d is the dipole matrix element 
of the resonant transition. |6L . The relative semplicity 
of the semi-classical derivation comes at the price of re- 
stricting the validity of the model to a well defined range 
of parameters: The concept of force is purely classical, 
therefore quantum fluctuations, stochastic heating and 
any incoherent process are to be neglected. This limits 
the validity of this model to large detunings | A| 3> w a , T 
(r is the atoms natural line width). Under these limi- 
tations, the potential Vd can be inserted into the atom 
Gross-Pitaevskii equation where it describes the laser- 
induced dipolc-dipole interaction between the atoms: 
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Here Ho is the linear single-particle Schrodinger Hamilto- 
nian, the wave function VP is normalized as N = / \^\ 2 dr 
with N denoting the total number of atoms, so that the 
gas density is n = \^f\ 2 , Uq = 4-7r?i 2 a s /m, m is the atom 
mass and a s is the s-wave scattering length (which will 



be assumed positive as for repulsive atom-atom interac- 
tions). Furthermore, since we are interested in the sta- 
tionary behaviour of the system and we have already as- 
sumed a stationary form for the electromagnetic field, we 
will consider ^(r,t) = $(r) exp(— iui a t). The atom equa- 
tion was already derived in [6[ within a fully quantum 
model and it is important to underline that, once the 
limitations of the semi-classical reasoning are taken into 
account, the two derivations lead to the same equation. 
To describe the role played by the electromagnetic field 
and the effect of the atoms on such field, it is necessary 
to include a field equation. Maxwell's equations for the 
propagation of radiation in a medium, 0, HH, [l4| ■> yield a 
wave equation which, under the assumption of L n 3> \l 
or Ve • E ~ (L n is the characteristic length scale of 
transverse density modulations and Al is the radiation 
wavelength), gives the three scalar equations (ojl = k^c) 
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When the input field distributions do no match the ex- 
act stationary solutions (which can be found numerically, 
0), propagation effects of some sort are to be expected. 
As demonstrated in 0, in the case of red detuning, the 
system settles down asymptotically to a stationary state 
with mutually localized atom-laser structures: Starting 
from a gaussian atom density profile and a super-gaussian 
laser intensity one, the interaction leads to the formation 
of two bell-shaped structures which propagate unchanged 
thereon. This means that atoms and radiation in excess 
will be shed away, which is the process we would like to 
elucidate here. Choosing z as the propagation coordinate 
and limiting the investigation to slow envelope variation, 
we consider 



£(r) = a(x, z) exp(ifc^2;)e, 
$(r) = ■4>(x,z)exp(ik a z), 



(3) 
(4) 



where x denotes the dimension transverse to the propa- 
gation direction z (one transverse dimension only for sim- 
plicity), e is the polarization vector of the field and k a is 
the atom wave number. The coupled system of equations 
(H|), ([21), can then be written in normalized variables as 
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where the following normalisation has been used: i = 
xkr,, for the atom wave function ip = ip/ip* with 
(47r|a|/3)?/' 2 = 1, for the laser a = a/a*, with 
ra|a|a 2 /(2?i 2 fc 2 ) = 1, s = signia), fi = k a /k^ (for 
simplicity we will assume fi = 1 hereafter) and (3 co u = 
6a s /(fc 2 |a|). The tilde will be dropped hereafter unless 
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otherwise stated. The red detuning case studied here will 
correspond to s = +1. Notice that no mutual localization 
is possible in the blue detuning case. While the classical 
description for the laser field is justified by the choice of 
the intensity regime, for a mean field model to be valid 
for the atom wave function, we must consider not only a 
zero temperature limit but also a low density limit with 
naj? <C 1, see Q. Furthermore, a low density regime is 
required in order to avoid the singularity of the model 
and consequent spurious collapse-like phenomena. 



III. INITIAL EVOLUTION 

As done previously, we will start assuming an initial 
laser intensity profile in super-gaussian form much wider 
than the gaussian initial atom density profile, both of 
them definitely different from the stationary solutions of 
the system thus ensuring a dynamical evolution: 
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where g is the supcrgaussian parameter (g = 10 in the 
simulations). The flat-top laser profile eliminates gradi- 
ent forces on the atoms at the very initial stage. How- 
ever, the flat top is immediately modified due to the nat- 
ural evolution of the system and the initial steps will be 
the seed of the subsequent structure generation. The 
atoms will imprint a chirp on the laser with the ef- 
fect of creating a central intensity peak with two lateral 
throughs, 0- This can be formally seen via a pertur- 
bative solution of the first propagation stage (i.e. for 
z <C With \%p\ 2 <C 1, the denominators in Eqs.© 

and ((6]) can be expanded keeping terms up to the order 
~ \ip\ 2 ■ Separating amplitude and phase as a(x,z) = 
A(x, z) exp(i6(x, z)),ip(x, z) = B(x, z)exp(i<fi(x, z)) the 
two equations give, upon separation of real and imagi- 
nary parts, 
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Consider a perturbative expansion F(x, z) 
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F 1 (x)z + F 2 {x)z 2 and G(x,z) = d(x)z + G 2 (x)z 2 up to 
second order in z where F stands for the functions A and 
B while G stands for 9 and <f> and the zero-th order terms 
are the initial functions (JT]) and (JSJ. Identifying powers 
of z, a solution is obtained for the amplitudes: 
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where /3'(x) = /3 — 2A 2 ) (x). This solution has the features 
observed in the initial evolution of the coupled system: 
The laser intensity profile changes in such a way as to 
peak in the center and at the same time two troughs 
are created on each side of the rising peak. The atom 
density profile shows the well known nonlinear defocus- 
ing behaviour - the center is depressed and two humps 
are created on both sides of the depression. This is the 
beginning of the creation of the stable mutually local- 
ized structures discussed in in a soliton-like process 
the nonlinearity in the atom equation can act as a self- 
generated trapping potential for the BEC. 



IV. STRUCTURE EMISSION 

As a consequence of the initial evolution stage, pro- 
vided the strenght of the focusing dipolc-dipolc interac- 
tion and that of the defocusing collisional nonlinearity are 
initially not completely out of balance, some atoms start 
to broaden away from the central structure while a part 
of the initial distribution remains trapped there. The ra- 
diation reacts to this process because of the dependence 
of the refactive index on the density profile and part of 
it is focused around the peak of the atom density. How- 
ever, if the trap induced by the laser is much wider than 
the atom wave function, the atoms lost from the central 
core can still be trapped. What initially was a hump of 
disperding atoms, can get trapped in a secondary self- 
generated potential well and induce mutual localization 
on the wings. The generation of these secondary mutu- 
ally localized structures, keeping the initial laser with and 
peak intensity fixed, should depend on having enough 
atoms escaping from the central peak since the escaping 
atoms must affect the laser wings to provoke the forma- 
tion of the secondary trap. Therefore we have numeri- 
cally studied the coupled evolution of (JTJl and ([5} for fixed 
d a = 5Ai,c?i = 8d a and fixed ao = 0.1346 correspond- 
ing to an initial peak laser intensity of 0.0153 mW/cm 2 , 
but varying i^q. In the simulations we have /3 ~ 38 cor- 
responding for instance to a detuning of 100 times the 
decay rate for 87 Rb atoms and s = +1. 
As anticipated, for low ipo, only a central density peak re- 
mains, a phenomenon studied in Q. The central atoms 
affect the laser profile which creates a trapping poten- 
tial. The atoms that escape this trap are not enough 
to modify the natural evolution of the laser wings which 
undergo well known modulations before diffracting away. 
Increasing ipo, the central structure generated by the sys- 
tem will be obviously modified, the balance of repulsive 
collisional interactions and attractive dipole forces has 
changed. Furthermore, the effect of escaping atoms be- 
comes stronger to the point that the same trapping mech- 
anism can now be realized on the sides of the central 
peak. FigO] shows the results of such an interaction for 
two values of ipo. The laser, as one would expect, forms 
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FIG. 1: Structure formation for (a) ipo = 0.0316 (correspond- 
ing to an initial peak atom density no = 1.7 10 19 m -3 ) and (b) 
ipo = 0.06645 (corresponding to no = 7.51 10 19 m~ 3 ). Dotted 
line: initial density distriubution. The propagation distance 
is indicated on the plots. All other parameters as specified in 
the text. All quantities normalized as in the text. 



analogous localized structures in correspondence of the 
atom density peaks. There is actually a formation of lo- 
calized structures even for low ipo ; the very low density of 
escaping atoms can focus extremely weak laser peaks, the 
process creates continuous families of mutually localized 
solutions. However, for low ipg they are hardly visible. 
What is interesting about these structures is their fate. 
They are self-consistcntly formed due to the effect they 
have on the laser radiation. Atoms focus the radiation, 
the radiation in turns exerts a focusing action on the 
atoms counterbalanced by their own dcfocusing interac- 
tion and their kinetic energy. During the initial transient, 
which lasts until the atom-laser structures are mutually 
adjusted to their own localized form, the lateral peaks 
are oscillating around the point where they have been 
trapped. They are kept there by the presence of the laser 
trap, laser wings have not yet completely adjusted to the 
newly born structures and they still act as an external 
trap for the atoms. Figs|2ja) and (b) show the inter- 
mediate stage of this transient for the same paramaters 
as in FigOJb). The structures are oscillating within the 
laser-induced trap which is being formed, FigsHJc) and 
(d), and once the laser has completely adjusted nothing 
keeps the atom-laser peaks oscillating around a fixed po- 
sition anymore and the structures are free to move away, 
Figs|2{e) and (f). For the parameters of Fig|2j they are 
ejected from the initial interaction region and proceed 
propagating with constant velocity as solitary-like waves. 
This could be explained by the repulsion due to the cen- 
tral peak: the two lateral peaks cannot proceed moving 
inward because they cannot overcome the repulsive bar- 
rier due to the central one. (???????) The position of the 
lateral peaks as a function of the propagation distance for 
the same parameters of Figj2]is shown in Fig|3|a), from 
which it is clear that, after an initial transient during 
which it is quite difficult to keep track of the structures' 
positions, the two "jets" are propagating at constant ve- 
locity. It is also evident how laser and atoms jets move to- 
gether. FigJ3{b) shows the peak value of the atom density 
of the emitted structures which tend to stabilize on a sta- 
tionary value. In a way, this phenomenon is reminiscent 
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FIG. 2: Details of the process of structure formation for 
ipo = 0.06645 (n = 7.51 10 19 m~ 3 ). Solid line: atom wave 
function, dotted line: laser-induced potential acting on the 
atoms (divided by 10 to make the figure more easily read- 
able). The propagation distance is indicated on the plots. 
All other parameters as specified in the text. All quantities 
normalized as in the text. 
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FIG. 3: (a) Emitted structure position versus propagation 
distance for ipo = 0.06645 (no = 7.51 10 19 m~ 3 ). Red: laser 
jet position; blue: atom jet position, (b) Peak atom density 
of the emitted structures for the same case. All quantities 
normalized as in the text. (Color on line.) 



of the emission of solitons engeneered in nonlinear op- 
tics with the aim for instance of implementing all-optical 
switching and directional couplers, (Toj . Whereas in the 
optics case the emission is stimulated only on one side, we 
obtain two moving structures because of the symmmetry 
of the configuration. We must underline that we refer to 
the emitted structures as solitary-like waves because of 
their ability to propagate with unchanged shaped but we 
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FIG. 4: Jet positions for different initial values of the atom 
peak density ipo- All other parameters are the same as for 
FigED (a) Vo = 0.052, (b) ip = 0.054, (c) Vo = 0.0662, (d) 
ipo = 0.0664, (e) ipo = 0.0668. All quantities normalized as in 
the text. (Color on line.) 



have not yet proved their collisional properties, prelimi- 
nary results indicate a behaviour strongly suggestive of 
a soliton-like nature. 

The analogy is also suggestive of the possibility of soli- 
ton steering. In fact, the properties of the structures 
ejected (peak density, velocity and number of jets) de- 
pend on the initial conditions. Therefore, changing the 
initial value of ipo , we have found jets emitted at differ- 
ent angles with respect to the propagation direction z 
and with different peak densities and peak laser intensi- 
ties, as can be seen from FigH] which shows jet positions 
for a few different cases. This last figure also shows the 
anomalous behaviour of the structures emitted starting 
from ipo = 0.0668. They initially move clearly inwards 
before being ejected. For growing initial peak density, 
there seems to be a stronger central trapping capable 
to attract the lateral peaks towards the center. Notice 
from FigQ] how, for higher initial ipo the jets tend to be 
born closer and closer to the central peak, where they 
are likely to experience a stronger interaction with it, 
due to a larger overlap (compare cases (a) and (b) in 
that figure with cases (c) and (d) which have larger ipo). 
There is a critical combination of parameters, which in 
our case occurs for ip = 0.0669, such that the two jets 
are drawn backwards until they collide and fuse at the 
center, FigJSJ It is known that the result of a collision 
between two solitons depending on the relative phase can 
lead to the fusion of the two objects, [l5| and references 
therein, however the nature of the collision within the 
model presented here needs further studies. After the 
merging, the remaining central peak stabilizes and does 
not undergo any dynamical changes anymore but it is 
very likely that such a structure will not be realized due 
to the extra-effects that are not considered within this 
model and that could play an important role during the 
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FIG. 5: Propagation and fusion for ipo = 0.0669 (no = 
7.62f0 19 m -3 ). Solid line: atom wave function, dotted line: 
laser-induced potential acting on the atoms (divided by 10 to 
make the figure more easily readable). Propagation distance 
as indicated on the plots. All quantities normalized as in the 
text. 



collision. For higher values of ipo no central peak is left 
while two lateral peaks are again symmetrically ejected. 
This could suggest an instability of the central peak as a 
possible explanation of the merging shown by the previ- 
ous case. If the central peak is unstable against diffrac- 
tion/defocusing and the laser-induced force is not able 
to keep it trapped, its atoms will broaden away with two 
possible outcomes for the jets: Either the repulsive in- 
teraction between the jets and the centrally disperding 
atoms is not strong enough to prevent the jets from merg- 
ing in the center, or it is important enough to push them 
away, compare Fig|5]and[5j 

It is interesting to notice how the integral of the jets 
wavefunction (N = \ip(x)\ 2 dx in FigJT]) seems to 
tend to a finite value as a function of the integral of the 
initial wavefunction (iVo in FigjT]). This would be ac- 
ceptable from the point of view of soliton behavior: The 
emitted solitary-like structures can accomodate a given 
number of atoms, atoms in excess will go and form extra 
jets, an example is shown in FigJS] 

A final note concerns one more analogy with an 
optical soliton behaviour. It seems in fact possible to 
excite a structure very similar to the bound system 
observed for optical solitons in which two pulses perform 
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FIG. 6: Structure emission for ?/>o = 0.092 (no = 
7.6210 19 m _3 ). Solid line: atom wave function, dotted line: 
laser-induced potential acting on the atoms (divided by 10 to 
make the figure more easily readable). Propagation distance 
as indicated on the plots. All quantities normalized as in the 
text. 
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FIG. 7: Integral of the jets wavefunctions N = \ip(x, z — > 

oo)| 2 dir, (solid line) and of the central peak wavefunction 
(dotted line) versus the integral of the initial wavefunction 
N — \ip(x,z = 0)| 2 da;. The points were the solid line 
is broken correspond to merging and fusion and therefore no 
emission of jets at all. All quantities normalized as in the 
text. 
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FIG. 8: Structures emitted for ip = 0.2 (n = 6.81 10 20 m" 3 ). 
Propagation distance as indicated on the plot. All quantities 
normalized as in the text. 




FIG. 9: Evolution of the central atom density calculated as 
\tp(x — 0)| 2 as a function of the propagation distance for ipo = 
0.196 (no = 6.54 10 20 m~ 3 ). All quantities normalized as in 
the text. 



an oscillatory motion b y b ouncing back and forth in 
their own potential well, [15[ • In a repeated dance, under 
particular conditions, the optical solitons pass through 
each other, move apart and come to a halt to move back 
together. This is what can be seen for a given choice of 
initial parameters for the system under analysis here. 
Figin] shows the value of the atom density at x = as 
a function of the propagation distance for ipo = 0.196 
and oscillations which would agree with the presence of 
a bound soliton state are quite evident. Corresponding 
snapshots are given in FigfTTJl 
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FIG. 10: Snapshots of the evolution of the atom wavefunction 
for ipo = 0.196 (no = 6.54 10 20 m -3 ). Propagation distance as 
given in the plots. All quantities normalized as in the text. 



V. CONCLUSIONS 

In conclusion, proceeding from the idea that laser-BEC 
dipole-dipole interactions can lead to mutually localized 
structures, we have analyzed in detail the mechanism of 
formation of such structures concentrating on the process 
through which the structures shed away the extra atoms 
and extra radiation. Numerical simulations seem to indi- 
cate the possibility of generating and emitting secondary 
solitary-like wave packets in a jet-like fashon. Although 
the model used here is strongly simplified and any com- 
parison with experiment will require major refinements, 
the equations we have used enlighten the main physical 
effects and it seems possible to choose parameter regimes 
in which the effects neglected here will not destroy these 
results. This processes could be a further evidence of 
the analogy between matter waves and optical waves and 
even open the discussion about applications such as soli- 
ton stirring in BECs. 
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